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Abstract 

The direct method of construction of local conservation laws of partial differential equations 
', (PDE) is a systematic method applicable to a wide class of PDE systems [Anco S. and Bluman 

G., Direct construction method for conservation laws of partial differential equations Part II: 
O ! General treatment. Europ. J. Appl. Math. 13, 567-585 (2002)]. Within the direct method, 

one seeks multipliers, such that the linear combination of PDEs of a given system with these 
multipliers yields a divergence expression. 

After local conservation law multipliers are found, one needs to reconstruct the fluxes of 
the conservation law. In this review paper, we discuss common methods of flux computation, 
compare them, and illustrate by examples. An implementation of these methods in symbolic 
software is also presented. 
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1 Introduction 

■ Local conservation laws of systems of partial differential equations (PDE), i.e., divergence expres- 

sions vanishing on solutions of a given PDE system, arise in a wide variety of applications and 
contexts. They often serve as mathematical expressions for fundamental physical principles, such 
as conservation of mass, momentum, charge, and energy. Conservation laws of time-dependent 
bounded systems describe conserved quantities. Conservation laws of ordinary differential equa- 
tions (ODE) yield first integrals. 
Conservation laws are widely used in PDE analysis, in particular, studies of existence, uniqueness 
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and stability of solutions of nonlinear PDEs (e.g., [1-3]). 

Another important area of application of conservation laws is the development and use of numeri- 
cal methods. In particular, knowledge of conservation laws of a PDE system of interest enables one 
to draw from a wide variety of robust and efficient numerical schemes which have been derived in 
the last decades, including finite volume, finite element, and discontinuous Galerkin methods (see, 
e.g., [4, Chapter II; 5, Chapter IV]). 

In addition, local conservation laws of a given PDE system are used to introduce nonlocal (po- 
tential) variables and thus construct PDE systems nonlocally related to a given one, with the same 
solution set. The related framework of nonlocally related systems (e.g., [6-11]) has yielded many 
new results over the recent years. 

When a PDE system has an infinite set of local conservation laws involving arbitrary functions, 
it can be sometimes mapped into a linear PDE system by an invertible transformation, which can 
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be derived explicitly when it exists [12]. An infinite countable set of local conservation laws is often 
associated with integrability. 

In the current paper, we consider local divergence- type conservation laws. Other types of con- 
servation laws (such as lower-degree, e.g., curl-type, and nonlocal conservation laws) also arise in 
applications. 

We now turn to the question of derivation of conservation laws. Many classical conservation laws 
have been found by ad hoc methods. One of the most well-known systematic methods is due to 
Emmy Noether, who demonstrated that for self-adjoint (variational) PDE systems, conservation 
laws arise from variational symmetries, i.e., symmetries that preserve the action integral [13]. 
However, the applicability of Noether theorem is severely limited, since the majority of PDE systems 
arising in applications are not self-adjoint [14]. 

An interesting review devoted to comparison of several methods of computation of conservation 
laws is found in [15]. 

A more general systematic method of constructing local conservation laws called the direct method 
was suggested in [16,17] (see also [18]). We review this method and the underlying theory in Section 
[3] of the current paper. Within this method, one seeks a set of local multipliers (also called factors 
or characteristics) depending on independent and dependent variables of a given PDE system and 
derivatives of dependent variables up to some fixed order, such that a linear combination of the 
PDEs of the system system taken with these multipliers yields a divergence expression. Families 
of multipliers that yield conservation laws are found from determining equations that follow from 
Euler differential operators. After finding sets of local conservation law multipliers, one needs to 
derive expressions for the corresponding conservation law fluxes. 

The goal of this paper is to review and illustrate available methods of flux construction (Section 
HI). We start from a direct ("brute force") method (Section 14. ip that converts a conservation law 
directly into the set of determining equations (linear PDEs) for the unknown fluxes. This method is 
easy to apply in the case of simple conservation laws, including those involving arbitrary functions. 

The second and the third methods (Sections 14.21 and 14. 3p employ integral formulas, related to 
homotopy operators, to compute the fluxes. These formulas can be used to produce closed-form 
flux expressions for rather complicated forms of conservation laws, except for the cases when the 
conservation law involves arbitrary functions. The first integral formula is presented following [19]. 
The second integral formula initially appeared in [17], but is re-derived in the present paper in a 
different form, which is symmetric with respect to independent variables, and is simpler for practical 
computations and software implementation. 

The fourth method [20] (Section 14.41) applies to scaling-homogeneous conservation laws of scaling- 
invariant PDE systems, which often occur in applications. Within this method, one finds fluxes of 
conservation laws through a rather simple formula involving no integration. 

Another appropriate reference for the above four methods is [21, Section 1.3]. 

In Section [5l we briefly compare and discuss the four flux computation methods presented in the 
current paper. 

Finally, in Section [6J we present a symbolic software implementation of the direct method of 
conservation law construction and the four flux computation methods in the symbolic package GeM 
for Maple [22], and illustrate these methods with examples. 

There also exist other symbolic software packages implementing the direct method, as well as 
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other methods for computation of conservation laws; see [19, 26-29] and references therein. 



2 Notation and definitions 

General definitions 

Consider a general system of N partial differential equations of order k with n independent variables 
X — i^X j ... j X ) and m dependent variables u(x) = (u (x), . . . ,u m (x)), denoted by 

R°[u] = R a (x,u,du,...,d k u) = 0, a = l,...,N. (2.1) 

Here partial derivatives are denoted by u? = 9u q^ ', 

du = d l u = (u\(x), ul(x), uf{x), it™ (a?) J 

denotes the set of all first-order partial derivatives; 

d p u = \v% iip | n=l,...,m; i x , . . . , i p = 1, . . . , n\ 

d p u^{x) 



dx 11 . . . dx l p 



fi = l,...,m;ii,...,i p = 1,... ,n 



denotes the set of all partial derivatives of order p. 

Within the current paper, the notation "/[it]" means that / is a function of one or more inde- 
pendent variables x, dependent variables it, and possibly derivatives of dependent variables, up to 
some fixed order, i.e., 

/[it] = f(x, u, du, . . . , d u), I > 0. 



Definition 1. The given PDE system (|2.ip is written in a solved form with respect to some leading 
derivatives, if 

R a [u]=u j ° aiias -G CT {x,u,du,...,d k u) = Q, a=l,...,N, (2.2) 

where s < k, 1 < j a < m, 1 < . . . , i a ,s < n for all a = 1, . . . , N. In (12. 2ft . {itf^ icr s } is a 
set of linearly independent sth order leading partial derivatives, with the property that none of 
them or their differential consequences appears in {C^w]}^^. (For example, for dynamical PDE 
systems, time derivatives yield a natural set of leading derivatives.) 



Definition 2. A PDE system (|2.ip is in Cauchy-Kovalevskaya form with respect to an independent 
variable x 3 , if the system is in the solved form for the highest derivative of each dependent variable 
with respect to x 3 , i.e., 

u a = G a (x,u,du,...,d k u), l<s {a) <k, a = l,...,m, (2.3) 



d{x3) 



where all derivatives with respect to x 3 appearing in the right-hand side of each PDE of (|2.3|) are 
of lower order than those appearing on the left-hand side. 
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Definition 3. A PDE system ()2.1|) admits a Cauchy-Kovalevskaya form if it can be written in 
Cauchy-Kovalevskaya form f|2.3|) with respect to some independent variable (after a point (contact) 
transformation if necessary). 

A PDE system can admit a Cauchy-Kovalevskaya form only if its number of dependent variables 
equals the number of PDEs in the system, i.e., N = m. Many PDE systems arising in applications 
admit a Cauchy-Kovalevskaya form with respect to one or more of its independent variables. PDE 
systems that do not admit a Cauchy-Kovalevskaya form are typically PDE systems with differential 
constraints, such as Maxwell's equations. 

Definition 4. The total derivative operators with respect to independent variables are given by 

D ^^ + <^ + <4 r+< -47 + ---' * =l, ---' n - (2 - 4) 



In (|2.4p and below, summation over repeated indices is assumed. 

Definition 5. The differential Euler operator with respect to a function U{x) is defined by 

Eu = ^--D~ + ... + f-im, . . . Di — — — + . . . . (2.5) 
u dU dUi v ; 1 s dU h ... is y J 

Local conservation laws 

Definition 6. A local conservation law of PDE system (|2.1[) is a divergence expression 

div$[u] =Di**[u] ^Di^u] + ... + D„,$ r >] =0 (2.6) 

holding for all solutions of PDE system (|2.1|) . In (|2.6p . = <I> J (x, it, <9u, . . . , cfit), i = 1, . . . , n 

are called the fluxes of the conservation law, and the highest order derivative (r) present in the 
fluxes is called the (differential) order of a conservation law. 



If one of the independent variables of the PDE system ([2.60 is time £, the conservation law ([2.60 
takes the form 

D t tf[u] + div*[u] = 0, (2.7) 

where div$[n] = Dj = Di$ 1 [n] + . . . + D n _i$ n_1 [u] is a spatial divergence, and x = 

(ar, . . . , x n ~ 1 ) are n — 1 spatial variables. Here ^[it] is referred to as a density, and as 
spatial fluxes of the conservation law ([2.70 . 



Definition 7. A local conservation law ([2.60 of the PDE system ([2.10 is trivial if its fluxes are of 
the form $*[«] = M*[u] + H x [u], where M l [u] vanishes on the solutions of the system ([2.10 . and 
T>iH l [u] = is a trivial divergence. 

Definition 8. Two local conservation laws ([2.10 are equivalent if they differ by a trivial conservation 
law. Similarly, a set of local conservation laws of the form ([2.10 is linearly dependent if its nontrivial 
linear combination (i.e., a linear combination with coefficients which are not all equal to zero) yields 
a trivial conservation law. 
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3 The direct method of construction of conservation laws 



For a given PDE system, one is interested in finding its nontrivial linearly independent local con- 
servation laws. 

We now outline a systematic and generally applicable method of construction of local conservation 
laws of PDE systems, called the direct method [16,17]. 

3.1 Multipliers. The sequence of the direct method 

Consider a set of multipliers {A a [U]}^ =1 = {A (T (x, U, dU, . . . , d l U)}^ =1 which, taken as factors in 
the linear combination of equations of the PDE system (|2.ip . yield a divergence expression 

K[UW\U) =Di$TO (3.1) 

which holds for arbitrary functions U{x). Then on the solutions U{x) = u{x) of the PDE system 
(|2.ip . one has a local conservation law 

A a [u]R a [u] = D;$>] = 0, (3.2) 

provided that each multiplier is non-singular. [A multiplier A CT [f7] is called singular if it is a 
singular function when computed on solutions U(x) = u(x) of a given PDE system (|2.1[) (e.g., 
A<j[i7] = F[U]/R a [U]). In practice, one is interested only in non-singular sets of multipliers, since 
singular multipliers can yield arbitrary divergence expressions which are not conservation laws of 
the given PDE system. For example, taking A CT1 [Z7] = div <&\U]/R ai [U] and A a ^ ai [U] = 0, one 
obtains A (T [i7]ii <T [l7] = div <&[£/], which does not yield a conservation law of the PDE system (|2.ip .] 

To seek sets of multipliers {A (T [C/]}^L 1 that yield conservation laws (|3.2p of a given PDE system 
(|2.ip , one uses the following fundamental property of Euler operators (|2.5p . 

Theorem 1. Let x = (x , ...,x n ) and U(x) = (U 1 (x), . . . , U m (x)). An expression f[U] = 
f(x, U, dU, . . . , d l U), I > 0, is annihilated by all Euler operators Eyj , i.e., 

E w f[U] = 0, j = l,...,m, (3.3) 
if and only if f[U] is a divergence expression 

f[U)=D i $ i {U], 
for some set of fluxes {&[U]}? =1 = U,dU, . . .,d r U)}f =l . 

Corollary 1. A set of non-singular multipliers {A (7 (x, U, dU, . . . , d l U)}^ =1 yields a local conserva- 
tion law for the PDE system (I2.ip if and only if 

E U3 (A a (x,U,dU,...,d l U)R a (x,U,dU,...,d k U)) = 0, j = l,...,m, (3.4) 

hold for arbitrary functions U{x). 
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The sequence of the direct method 

1. For a given PDE system (|2,ip . seek multipliers of the form {A a (x, U, dU, . . . ,d l U)}^ =1 , for 
some specified order I. Choose the dependence of multipliers on their arguments so that 
singular multipliers do not arise. 

2. Solve the set of determining equations (|3.4p for arbitrary U(x) to find all such sets of multi- 
pliers. 

3. Find the corresponding fluxes $ l [U] = & l (x, U, dU, . . . , d r U) satisfying the identity (|3.1|) . 

4. Each set of fluxes and multipliers yields a local conservation law 



holding for all solutions u(x) of the given PDE system (|2.ip . 

In practical computations, when it is possible, it is convenient to write the given PDE system (12. ip 
in a solved form (|2.2p with respect to some leading derivatives. Then in Step 1 of the direct method, 
one can easily avoid singular conservation law multipliers, by excluding such leading derivatives 
and their differential consequences from the multiplier dependence. 

3.2 Completeness of the direct method 

The direct method outlined above can be applied to a wide class of PDE systems. Indeed, the only 
requirement is the existence of all derivatives that arise in the determining equations (|3.4[) . 

However, for a given PDE system, it is often important to know whether or not the direct method 
can yield all its nontrivial local conservation laws (at least up to some fixed order). Conversely, 
under what conditions does a set of local multipliers yield only nontrivial local conservation laws? 
We now outline the results related to these questions. 

When a given PDE system is written in a solved form, one can show that the application of the 
direct method with multipliers of sufficiently general dependence yields all local conservation laws 
of the given PDE system of any fixed order order I. In particular, the following theorem holds [21]. 

Theorem 2. For each local conservation law Dj<J> l [ii] = of a PDE system (I2.ip written in a solved 
form ()2.2p . there exists an equivalent local conservation law Dj<E» l [u] = which can be expressed in 
a characteristic form 



in terms of a set of non-singular local multipliers {A (T [?7]}^L 1 , with fluxes that contain no leading 
derivatives U-^ ia g nor their differential consequences. 

However, even when a PDE system is written in a solved form (12. 2p . the correspondence between 
equivalence classes of conservation laws and sets of conservation law multipliers might not be one- 
to-one. The following theorem holds [17]. 



D,$ l (x,u,<9u,...,d r u) =0, 




(3.5) 
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Theorem 3. Suppose a PDE system (|2.ip admits a Cauchy-Kovalevskaya form (|2.3p . Then all 
of its non-trivial (up to equivalence) local conservation laws arise from multipliers. Moreover, 
there is a one-to-one correspondence between equivalence classes of conservation laws and sets of 
conservation law multipliers {A a [U]}^ =1 , if the multipliers do not involve derivatives with respect 
to xK 

4 Computation of fluxes of conservation laws 

Consider a set of non-singular multipliers {A CT [{7] = A CT (x, U, dU, . . . , d l U)}^ =1 that yields a diver- 
gence expression (|3.ip for a PDE system (|2.ip for some set of fluxes = U, dU, d r U)}f =1 
Hence one has a local conservation law (|3.2|) holding on the solutions U(x) = u(x) of the PDE sys- 
tem flzrp . 

Given a set of conservation law multipliers, the problem of finding the fluxes {<]?*[£/] }™ =1 is formally 
a problem of inversion of the divergence differential operator. Even modulo trivial conservation 
laws of the first type (Definition H l [u] = 0) , a set of multipliers defines an equivalence class 
of conservation laws, up to free constants and arbitrary functions. For example, suppose in a 
three-dimensional space (x,y,z) (n = 3), the fluxes of a conservation law 

X;A (7 [C/] j R ct [C/]^^D^[C7] 

cr=l i=l 

are given by 

$[U] = (<!> 1 [U],&[U\,$ 3 W) = 6iv- 1 ^J2A a [U]R a [U]J . 

Then the expression 

*[£/] = *[[/] + curl K, 

for an arbitrary smooth vector field K, also yields fluxes of an equivalent conservation law. From 
the practical point of view, for a given set of multipliers, it is sufficient to find any one set of 
corresponding fluxes. 

We now consider different ways of finding fluxes of conservation laws from known multipliers. 
4.1 The direct method of flux computation 

For conservation laws arising from simple forms of multipliers, the fluxes are most easily found by 
direct matching of the two sides of equation (|3.1|) and subsequent integration by parts. 

If integration by parts is not obvious, one can assume a sufficiently general dependence for the 
fluxes and look for a particular solution of PDEs (|3.1|) for the fluxes within such an ansatz. 

Here it is important to note that if the maximal order of derivatives of U(x) present in the multipliers 
is I, and the maximal order in the equations R a [U] appearing in the linear combination (|3.ip is 
k, then, without loss of generality (i.e., up to trivial conservation laws) one may assume that 
the fluxes are of the form <J? l [[7] = U, dU, . . . , d r U), where r = max(Z,fc) [23]. Using the 
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independence of & l [U] from higher derivatives, one can subsequently split PDEs (|3.ip into a set of 
simpler determining equations. 
We now consider an example of a direct flux computation. Further examples are found in Section 

El 

Example. As an example, consider the nonlinear wave equation 

R[u} = u tt -{c 2 {u)u x ) x = 0, (4.1) 

for arbitrary wave speed c[u). For simplicity, consider multipliers of the form A[U] = A(x,t,U). 
The determining equations f|3.3[) yield the solution A(x, t) = Ci + C^x+Cai+C^xi, where C±, . . . , C4 
are arbitrary constants. As a result, one obtains four linearly independent conservation laws, arising 
from the multipliers A" = 1, A^ 2 ) = x, A^ = t, A*- 4 ) = xt. We now determine the corresponding 
density- flux pairs. 
For the multiplier A^ 1 ) = 1, one obviously has 

AW[U]R[U] = B t (U t )-B x (c 2 (U)U x ), (4.2) 

since PDE (|4.ip is in divergence form as it stands: 

D t (u t )-D x (c 2 (u)u x ) = 0. (4.3) 

For the multiplier A^ 2 ) = x, one can determine the flux and density using integration by parts: 



AW[U]R[U] = x(D t (U t )-B x (c 2 (U)U x )) 

= D t (xU t ) - D x (xc 2 (U)U x ) + c 2 (U)a 



B t (xU t )-B x (x<?(U)U x -J . 



(4.4) 



c 2 {U)dUj , 

and thus the corresponding conservation law is given by 

B t (xu t ) - B x (^xc 2 (u)u x - j c 2 (ii)rfu^ = 0. (4.5) 

Similarly, for the multiplier A® = t, one finds the corresponding conservation law 

B t (tut - u) - D x {tc 2 (u)u x ) = 0. (4.6) 

To find the flux and density of the somewhat more complicated fourth conservation law of PDE 
(|4.1|) following from the multiplier A^ 4 ) = xt, one needs to solve the flux-density determining 
equation 

AW[U]R[U] = xt{D t (U t )-D x (c 2 (U)U x )) = D t ^[U]+D x ^[U}. (4.7) 

Since the left-hand side of (14. 7|) is linear in the highest derivatives Uu and U xx , one can assume 
that *[[/] = V(x,t,U,U x ,U t ) and $[£/] = $(x,t,U,U x ,U t ). Expanding both sides of lf£7]> . one 
obtains 

xt{U tt - 2c(U)c'{U){U x ) 2 - c 2 {U)U xx ) 

= {% + ^uU t + V Ut Utt + VuJJxt) (4-8) 
+ ($ a + <$>uU x + 3>u t U xt + $u x U xx ) . 
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Matching the terms of the highest order derivatives Uu, U xx and U x t, one finds that 

xt = ¥ Ut , -xtc 2 (U) =$ Ux , * Ux = -*u t . (4.9) 

The third equation in (|4.9p can be replaced by ^u x = <&u t = 0> since other choices lead to equivalent 
conservation laws. The first two equations in (14. 9|) yield 

V[U]=xtU t + a(x,t,U), $[U] = -xtc 2 {U)U x + p{x,t,U), (4.10) 

for arbitrary a(x,t,U), ft(x,t,U). Substituting (|4.1U|) into the determining equations (|4.8|) and 
setting to zero coefficients of first-order partial derivatives of U, one finds 

x = -au, tc 2 (U) = Pu, a t = -(3 X . 

Therefore 

a(x,t,U) = -xU + a(x,t), f3(x,t,U) = t j c 2 (U)dU + 0(x,t), a t = -/3 x . 

It is evident that any choice of a and (3 satisfying at = —(3 X yields an equivalent conservation law, 
with the simplest one having a = (5 = 0. Thus, the fourth conservation law of PDE (|4.ip is given 
by 

D t (xtut - xu) - D x ^xtc 2 (u)u x - t J c 2 (u)di?j = 0. (4.11) 
4.2 The first homotopy formula 

In the case of complicated forms of multipliers and/or equations, for the inversion of divergence 
operators, one can use homotopy operators that arise in differential geometry and reduce the 
problem of finding fluxes to a problem of integration in single-variable calculus. In this section, we 
present the first set of such formulas, following [19]. 

Definition 9. The n- dimensional higher Euler operator with respect to a function U(x 1 , ...,x n ) is 
given by 

E & *' ) = t - t ( M - ( K ) °"r" . . ■ »t--- n ° . (4.i2) 



where one denotes 

Qkl+.,.+kn,JJ 



si 



jj(ki+...+k n ) 



X" 



Note that E^' ^ = Efj is the usual Euler operator (|2.5p with respect to U. 

Definition 10. Let x = (x 1 ,...,x n ) and U = (U 1 (x), . . . , U m (x)). For an expression f[U] = 
f(x,U,dU,...), the n- dimensional homotopy operator is defined through its n components corre- 
sponding to n respective independent variables x l , i = 1, . . . , n, as follows: 

R^(f[U])= / Y.^VftVf, (413) 
Jo ~[ u=xu A 
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with 

f\fp]) = E-E + J , , 

Sl=0 s„=o Vi + Sl + '" + Sn/ (4.14) 
xDf ...D^ ^iE^-'* 1 - 1 ' 1+84, " +1 -" , ' B) (/[t?])) , 

j = l,...,m. 

The following theorem holds [19]. 
Theorem 4. Suppose f[U] is a divergence expression 

f[U] = div&[U] = D^p] + . . . + D n $ n [U], 
and f[0] =0. Then the fluxes $> X [U] are given by 

&[U]=R&)(f[U]), i = l,...,n. (4.15) 
up to corresponding fluxes of a trivial conservation law, provided that the integrals (|4.13p converge. 

Two examples of the use of the first homotopy formula are now considered. Further examples are 
found in Section [6l 

Example 1. As a first example, consider the Korteweg-de Vries (KdV) equation 

u t + uu x + u xxx = 0. (4.16) 

One can show that PDE (|4.16|) has a local conservation law arising from the multiplier A[U] = U, 
i.e., 

f[U] = U{U t + UU X + U xxx ) = DMU] + D X $[U]. 

is a divergence expression. We now reconstruct the corresponding density ^[U] and flux <J>[?7]. 
First compute 

= D t °D° (uE^°\f[U])) = d\ 



I^(f[U}) = D°D 2 (uEf 3 \f[U])) + D°Di (uEp 2) (f[U])) + D°D° (uE^ (f[U]) 
U 3 - Ul + 2UU XX . 



Consequently, 



H<*>(/[E7]) = (A 3 C/ 3 - A 2 [/ 2 + 2A 2 C/^) ^ = if/ 3 - it/J + LTCU 
J o A 

The corresponding conservation law is given by 

D t (|n 2 ) + D x (±w 3 - \u\ + uu«) = 0. (4.17) 
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Example 2. Consider the 2-dimensional G-equation 



R[g] =9t~\ grad£f| = g t - ^ g 2 x + g 2 = 0, (4.18) 

which describes flame propagation in a static gas (see, e.g., [24]). Here g(t,x,y) = implicitly 
defines the position of the flame surface at time t, and the surface advances at a constant speed in 
the normal direction. Let (x l ,x 2 ) = (x,y). In particular, one can show that the PDE (|4.18p has 
the following multipliers of conservation laws: 

A«[G] = -4 (G x G yy - G y G xy ), A®[G\ = H(G x ,G y )(G xx G yy - G%), (4.19) 

(j y 

where H(G x ,G y ) is an arbitrary function. For more details, see [25]. 
First consider the multiplier A^^G] in (|4.19p . From (|4.14p . one obtains 

_G 



I^(A« [G]R[G]) = ^(G x G xy - G y G xx ) 



which is an invariant expression under the transformation U = XU. The same is true about 
expressions (A^ [G]R[G]) and 1^ (A^ [G]R[G]) . Hence the integral homotopy operators ([4TT3"]) 
yield a divergent integral dX/X. 
Second, consider the multiplier A^ 2 ^[G] in (|4.19p . Since the corresponding divergence expression 



f[G] = A®[G]R[G] = H(G x ,G y )(G t - JG 2 X + G 2 y )(G xx G yy - G 



; 2 ^ 

T X y ) 



contains an arbitrary function H(G x ,G y ), the corresponding integrals (|4,13p cannot yield closed- 
form expressions. 

It follows that the present flux construction method cannot be used for either of the multipliers 
dUSD of the G-equation (gUJ. 

4.3 The second homotopy formula 

An alternative set of integral formulas for fluxes of local conservation laws was developed by Bluman 
and Anco [17]. Here we present these formulas in a different form, which is symmetric with respect to 
independent variables, and also is simpler for practical computations and software implementation. 

It takes a somewhat bigger effort to implement the current method in symbolic computation 
software than to implement the method presented in Section 14.21 However the formulas given in 
the current section have two important practical advantages over the ones presented in Section T4.2i 



(i) They can be applied to divergence expressions f[U] with /[0] ^ 0. 

(ii) They have flexibility that can be used to always avoid divergent integrals. 

The following theorem holds. 

Theorem 5. Suppose a PDE system (I2.ip has a local conservation law arising from a set of 
multipliers {A cr [f7]}^ r =1 , i.e., ()3.2p holds. Then the fluxes {§ l \U]}f =l of that conservation law are 
given by 

&[U]= ¥[U\+ f (S i \u-U, A[XU + (1 - X)U}; R[XU + (I - X)U] 

Jo K 1 J (4.20) 

U-U, R[XU + (1 - X)U];A[XU + (1 - X)U] ) dX, i = 1, . . . , n, 
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up to fluxes of a trivial conservation law. In (|4.20p . operators S l , S l are defined by their action on 
arbitrary functions V = {V 1 {x) , . . . ,V m {x)) , W = {W x {x), . . . , W N (x)), V = {V l (x), ... , V m (x)) 
and W = (Wi(x), . . . , Wn(x)) as follows: 

S*[V,W;R[U]]=Y, E (Dfr ■ ■ ■ D^W) D A . . . D jq W 9 ™ [UJ ) , (4.21) 

p=0 q=0 V J . •••./■/":••• 

A[E7]] = E E K • • • D vH % • • • f^ % r/ p A{r[ ^ ) . ( 4 - 22 ) 

p=0 g=0 V h-3qih- 

where k is the order of the given PDE system R{x ; u} (|2.ip . Z is the maximal order of the derivatives 
appearing in the multipliers [A^lU]}^^, and j\ . . . j q , i± ... i p are ordered combinations of indices 
such that 1 < jx < . . . < j q < i < i\ < . . . < i p < n. 
In (|4.20p . U = U(x) is an arbitrary fixed function, and & l [U] are quantities satisfying 

= A a [U]R a [U] =: F{x). (4.23) 

The proof of Theorem [5] is presented in Appendix EJ 

In ()4.20p . U is an arbitrary function of x, chosen so that the integral (|4.20p converges. Different 
choices of U yield fluxes of equivalent conservation laws. One normally chooses U = (provided 
that the integral converges). Once U = U{x) has been chosen, the functions {^[l/]}^ can be 
found by direct integration from the divergence relation (|4.23p by integrating F(x). For example, 
one can choose 



= J F(x)dx 1 , $ 2 [U] = . . . = $ n [U] = 0. 



It is important that unlike the case with the first homotopy formula (Section I4.2p . the proof of 
Theorem [5] does not require that the divergence expression A cr [C/]i2 CT [Z7] vanishes for U = 0. 

The formula (|4.20p is useful in many practical situations when the integrals can be explicitly 
evaluated. 

Two examples of the use of the second homotopy formula are now considered; further examples 
are found in Section [6j 

Example 1. As a first example, consider again the KdV equation (14.161) with multiplier A[U] = U. 
Using [7 = 0, one finds $ 1 [0] = <I> 2 [0] = 0, and hence the formulas (I4.20P yield exactly the same 
fluxes as the first homotopy method (Section 14. 2p . i.e., the conservation law (|4,17p . 
Alternatively, using U = x, one finds ( I )1 [C/] = tx 2 , $ 2 [C7] = 0, and hence 

<$> l [U]=tx 2 + \(U 2 -x 2 ), 

&[U] = I - Ix 3 + if/ 3 - \V 2 X + UU XX , 
which are a density and a flux of a conservation law equivalent to (|4.17p . 

Example 2. Consider the 2-dimensional G-equation (|4.18p and its local conservation law mul- 
tipliers (|4.19p . For the multiplier AW[Gr], it again yields a divergent integral when U = 0. The 
singularity can be removed by taking e.g. U = x, however even this simple choice leads to an 
integral with a highly complicated integrand. For the multiplier A^ 2 ) [G] involving an arbitrary 
function, formulas (|4.20|) do not yield a closed-form expression. 
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4.4 Computation of fluxes of conservation laws of scaling-invariant PDE sys- 
tems 



Consider a PDE system that has a scaling symmetry 

X s M= p (VA + ,(P) n pJL (4.24) 

where p^ and are constant scaling weights of independent and dependent variables, respectively, 
i.e., the generator (|4,24p corresponds to the global transformation group 

x i _> = e £X 4«l x * — e ep(l) x* 

u p _» u p = e £X ^u p = e eq[p) u p . 

The scaling symmetry generator (I4.24p can be written in evolutionary form as 

*.[«] = V[u]^ p = (>V -p<V<) A. (4.25) 

Assuming that the given PDE system can be written in a solved form (|2.2p . we first note that the 
scaling homogeneity of each PDE of the given system follows from considering the action of the 
scaling symmetry (14.24p on the leading derivative in each R a [u]. In particular, one has 

xW[U]R a [U] = r^R a [U], 

where A a ' = const is the scaling weight of each PDE R a [U]. 
Now suppose the given PDE system (12. 2p has a conservation law 



K[u]R a [u] = D f $*[u] = 0, (4.26) 

with fluxes = &{x,U,dU, . . . ,d l U). Suppose that this conservation law is scaling-invariant, 

and moreover, homogeneous under the scaling symmetry (|4.24|) . i.e., 

xW[J7]Di = PBi &[U], (4.27) 

where Xg is a corresponding prolongation of the scaling symmetry generator (|4.24|) , and P = const 
is a scaling weight of the conservation law. 

Then one can show [21] that each of the multipliers Ao-[{7] appearing in (|4.26p is homogeneous 
under the scaling symmetry (|4.24j) . In particular, X^[i7]A CT [[/] = sr a )A a [U], where 8r a \ = P — 
is the scaling weight of each Ao-[C/] ^ 0. 

Definition 11. A conservation law (|4.26p invariant with respect to a scaling symmetry X s [«] (|4.24p 
is called noncritical with respect to X s [tt] if 

n 

i=i 

for each a such that A CT [{7] ^ 0. 
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The formula for fluxes 

It follows that fluxes of noncritical homogeneously scaling conservation laws of scaling-invariant 
PDE systems can be found through a simple formula that involves no integration! The following 
theorem holds [20,21]. 

Theorem 6. Suppose the PDE system in the solved form (|2.2p has a scaling symmetry (|4.25|) . 
and a conservation law (|2.6|) . which is homogeneous under the scaling symmetry (|4.24|) . and is 
noncritical with respect to it. Then the fluxes & l [U] of such conservation law are given by 

**M=£ E (~l) g (Dfr ■ ■ ■ T> ip V) Da ■ ■ ■ D ig ^ [17] M , (4.29) 

p=0 q=0 V 31—jqm—tp/ 

up to fluxes of a trivial conservation law. In (|4.29p . is i/ie maximal order of derivatives appear- 
ing in the PDE system (| 2 . 1 [) . / is i/ie maximal order of derivatives appearing in the multipliers 
{A a [U]}^ =1 , and j\ . . . j q and i\ . . .i p are ordered combinations of indices such that 1 < j\ < . . . < 
jq < * < ii < • ■ ■ < < n. 

Remark 1. From the proof of Theorem [6] it follows that the formulas (|4.29p does not yield fluxes of 
critical conservation laws. In particular, for critical conservation laws, formulas (14.290 yield fluxes 
of trivial conservation laws. 

In practice, if a PDE system has two or more scaling symmetries, it can be the case that a 
particular conservation law is critical with respect to one scaling and noncritical with respect to 
the other, as is the case in the Example 2 below. 

Remark 2. For first- and second-order PDE systems, the flux formula (|4.29p takes a particularly 
simple form [20]. For second-order PDE systems, it becomes 

&[U] = fjP[U}A a [U}^^ + (B s nU])AAU]^P - V p Pn ( A ^]^P) > ( 43 °) 

where the summation is taken over repeated indices satisfying s > i and j < i. For first-order PDE 
systems, it further reduces to 

^[U]=rj p [U]Aa[U]^P- (4-31) 



Remark 3. In general, one can obtain fluxes of conservation laws using formula (|4.29|) not only 
for scaling symmetries, but for any pair consisting of a local symmetry X s [n] = i?[^]gfp (in the 
evolutionary form) and a set of local conservation law multipliers {A (T [C/]}^L 1 . Details are given in 
Appendix iBl (see also [16,21]). 

Example 1. As a first example, we consider again the KdV equation (|4.16p with multiplier 
A[U] = U. The PDE (|4.16j) has a scaling symmetry 
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Using the formula (|4.29p with {x ,x ) = (t, x), one obtains the density and the flux 

^[U] = {2U + 3tU t + xU x )U, 

&[U\ = U{3tU txx + xU xxx ) - 3tU x U tx 

+3{2U + tU t )U x x + U 2 (2U + 3tU t ) + U x (xU 2 - 3U X ), 

i.e., a conservation law 

B t ( (2u + 3tu t + xu x )u] + B x ( u(3tu txx + xu D3; ) - 3tu x u tx 

v y v v (4.33) 

+3(2u + tut)u x x + u 2 {2u + 3tut) + u x (xu 2 — 3u x ) J = 0, 

equivalent to the conservation law (|4,17p . 

Example 2. Consider the 2-dimensional G-equation (|4.18p and its local conservation law multi- 
pliers (|4.19|) . The G-equation (|4.18|) is a first-order PDE with obvious scaling symmetries 

d d d d 

X 1 = 5 — , X 2 = t— + x— + y— . (4.34) 
og ot ox oy 

Denoting (x 1 ,x 2 ,x 3 ) = (t,x,y), for Xi, one has q = 1, = = p^ = 0; for X 2 , q = 0, 
p (l) = p (2) _ p (3) _ ! 

Consider the conservation of the G-equation (|4.18p law with multiplier A^[G] = -^{G x G yy — 
G y G xy ) (I4TT91) . The function A^fG] is evidently homo geneous under both symmetries Xi and X 2 
(|4.34p . with respective scaling weights s% = — 1 and s 2 = 0. The corresponding scaling weights of 
the G-equation (|4.18p are r± = 1 and r 2 = —1. 

Using Definition 1111 one finds \i = ~ 1 + 1 = 0, and X2 = 0— 1 + 3 = 2, hence the conservation 
law with multiplier [G] is critical with respect to the scaling symmetry Xi and noncritical with 
respect to X 2 . Indeed, one can directly check that the conservation law obtained by formulas (|4.29p 
using the symmetry Xi, 



gl+alj \\9x + 9y 



0. 



is trivial. The correct conservation law fluxes corresponding to the multiplier A^ 1 ) [G] can be found 
using the symmetry X 2 ; such conservation law is given by 



fj[g] = -tgt - xg x - yg y . 
5 Comparison of the four flux computation methods 

In Table [lj we briefly outline the characteristics of the four methods of flux computation discussed 
above. 
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Table 1: Comparison of Four Methods of Flux Computation of Section [U 



Method 


Applicability 


Computational complexity 


Direct 

(Section HU) 


Simpler multipliers/PDE systems, which may in- 
volve arbitrary functions. 


Solution of an overdeter- 
mincd linear PDE system 
for fluxes. 


Homotopy 1 
(Section l4~2f 


Complicated multipliers/PDEs, not involving arbi- 
trary functions. 

The divergence expression must vanish for U = 0. 
For some conservation laws, this method can yield 
divergent integrals. 


One-dimensional integra- 
tion. 


Homotopy 2 
(Section P) 


Complicated multipliers/PDEs, not involving arbi- 
trary functions. 


One-dimensional integra- 
tion. 


Scaling sym- 
metry 

(Section OJ 


Complicated multipliers/PDEs, may involve arbi- 
trary functions. 

Scaling-homogeneous PDEs and multipliers. 
Noncritical conservation laws. 


Repeated differentiation. 



One can see that there is no "preferred" method of finding fluxes that is simple to use and 
sufficiently general to be recommended for a generic conservation law of a nonlinear PDE system. 
The following observations and recommendations are based on the author's experience. 



• If a given PDE system and multipliers do not involve arbitrary functions, one usually uses 
one of the two integral methods. Though the second integral method is more general, the 
first one is also sometimes used because it can yield simpler flux expressions. 

• If a given PDE system and/or multipliers are not very complicated but involve arbitrary 
constitutive functions, one normally attempts the direct method of flux computation. 

• When a given PDE system is scaling-invariant and the conservation law is scaling-homogeneous, 
then the scaling symmetry method is the method of choice, since it involves simplest com- 
putations. (However, as seen in the KdV example (|4.33p . when some ^ 0, the densities 
and fluxes computed through the fourth method can contain multiple unnecessary terms 
corresponding to trivial conservation laws.) 

• Another instance when one needs to use the fourth method is the case of complicated scaling- 
homogeneous PDEs and/or multipliers involving arbitrary functions. In this situation, the 
scaling symmetry method is the only systematic method available. 

In the following Section [5j we consider a software implementation of the four above-discussed flux 
computation methods, applied to several practical examples. 
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6 Symbolic software implementation of flux computation meth- 
ods. Examples 

The direct method of construction of conservation laws (Section 13. ip is algorithmic and has been 
implemented in symbolic softwarj^l. Here we discuss its Maple implementation in the symbolic 
package GeM [22] written by the author. The package GeM also contains routines for the four 
methods of flux reconstruction considered in this paper. 

6.1 General conservation laws of the nonlinear wave equation 

Here we seek local conservation laws for the nonlinear wave equation (|4.ip arising from the multi- 
pliers of the form A[U] = A(x,t, U) for an arbitrary c{u). 
One defines the variables and the PDE (|4.ip in GeM using the following Maple commands: 

with(GeM) : 

gem_decl_vars(indeps=[x,t] , deps= [U(x,t)] , f reefunc= [C(U(x,t) ) ) ; 
gem_decl_eqs( [diff (U(x,t) ,t,t)=diff (C(U(x,t) ) ~2*dif f (U(x,t) ,x) ,x)] , ' *' 
solve_f or=[diff (U(x,t) ,t,t)] ) ; 

The option solve_f or is used to specify a set of leading derivatives the given PDE systems can 
be solved for. Note that the direct method of conservation law computation does not require the 
equations to be written in the solved form, and therefore the specification of the option solve_f or is 
not necessary. However in (16. If) . this option solve_f or was specified, since it will be used later in the 
flux computation routine, which automatically verifies the correctness of flux computations, which 
is accomplished by explicitly checking that the conservation law divergence expression vanishes on 
solutions U(x) = u(x) of the given PDE (|4.ip . 

The set of determining equations for the local conservation law multipliers is obtained and sim- 
plified using the routines 

det_eqs : =gem_conslaw_det_eqs( [x,t , U(x,t)] ) : 
CL_multipliers : =gem_conslaw_multipliers () ; 

(6.2) 

simplif ied_eqs : =DEtools [rif simp] (det_eqs , 

CL_multipliers , mindim=l, arbitrary= [C(U)] ) ; 

The first command in (|6.2p sets up the set of local conservation law multipliers determining equa- 
tions (|3.4p . and splits them. The splitting is done using the fact that the determining equations 
(|3.4p are polynomial expressions in terms of derivatives U of orders four and higher, and such 
derivatives are linearly independent. This yields an overdetermined system of 7 determining equa- 
tions for the unknown local multiplier A[t/]. After using the rif reduction algorithm, the system 
reduces to three PDEs 

d 2 A[U] _ d 2 A[U] _ dA[U] 
dx 2 ~ dt 2 dU ' 
The four obvious linearly independent solutions of these determining equations can be obtained 
using the command 

multipliers_sol : =pdsolve (simplif ied_eqs [Solved] ) ; 



a ^In particular, in the software by T. Wolf [26] written for REDUCE computer algebra system. 
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and are given by 



A[U] =d + C 2 x + C 3 t + C 4 tx 



where C\, . . . , C4 are arbitrary constants. 
Finally, we apply the direct method of flux computation (Section I4.ip . using the command 

gem_get_CL_f luxes (Lam_sol , method="Direct") ; 

which generates the divergence expressions of the four conservation laws (14. 3D . (|4.5f) . (|4.6f) and 



6.2 The first four conservation laws of the KdV equation 

As it is well-known, the KdV equation (|4.16p has an infinite countable sequence of local conservation 
laws of increasing orders. Here we compute the first four of these conservation laws, by limiting 
the multiplier dependence, as follows: 



The variables and the PDE (|4.16p are defined using the commands 
with(GeM) : 

gem_decl_vars(indeps=[x,t] , deps= [U(x,t)] ) ; 

(6.3) 

gem_decl_eqs( [diff (U(x,t) ,t)=U(x,t)*diff (U(x,t) ,x)+dif f (U(x,t) ,x,x,x)] , 
solve_for= [diff (U(x,t) ,t)] ) ; 

The set of determining equations for the multipliers is obtained and simplified using the routines 

det_eqs : =gem_conslaw_det_eqs( [x,t , U(x,t), dif f (U(x,t) ,x) , 
diff (U(x,t) ,x,x)]) : 

(6.4) 

CL_multipliers : =gem_conslaw_multipliers () ; 

simplif ied_eqs : =DEtools [rif simp] (det_eqs, CL_multipliers , mindim=l) ; 
and is solved using the Maple command 

multipliers_sol:=pdsolve (simplif ied_eqs [Solved] ) ; 

to yield the four local multipliers 



A[U]=A(t,x,U,U x ,U xx ). 




(6.5) 



where Ci, . . . , C4 are arbitrary constants. 
We now reconstruct the fluxes. 



1. The first homotopy formula. Using the command 



gem_get_CL_f luxes (Lam_sol , method="Homotopyl") ; 



(6.6) 
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we call the routine for the first homotopy method described in Section B~2l It yields the conservation 
laws 

Df(n) + D x .(±u 2 + u xx ) = 0, 
l„,2\ , n I l„,3 1„,2 



D t Uu'\ +DJi!« d - ±u x + uu xx \ = 0, 

/ \ / , , (6-7) 

D t ( xu - \tu 2 \ + D X ( - \tu 3 + i(xu 2 + tu 2 x ) - u x + (x - tu)u xx J = 0, 

D t ( lu 3 + \UU XX ) + DJ ±U 4 + - UU tx + ^Waaj + u xx y\ = 0. 



2. The second homotopy formula. The corresponding routine for the method described in 
Section 14.31 is called as follows 



gem_get_CL_f luxes (Lam_sol , method="Homotopy2") ; 

and yields the four divergence expressions in the same form as in (|6.7|) . By default, the arbitrary 
function U = is used. One can override this setting and use, e.g., U = x, by using an option, as 
follows: 

gem_get_CL_f luxes (Lam_sol , method="Homotopy2" , arb_f unc_Homotopy2={U=x}) ; 



which yield four conservation laws equivalent to the ones given in (|6.7p . but having somewhat more 
complicated fluxes: 

D t (x(t - 1) + u) + D x (±(u 2 - x 2 ) + u xx ) = 0, 

2 ± l/„,2 ^2\\ i ri I I 1^.3 i 1„,3 1„,2 

xx I — 



D t [tx 2 + ±(u 2 - x 2 )J + D X ^- ±x 6 + ±u 6 - \u x + uu 
D t ( \x 2 (3t - t 2 - 2) + xu - \tu 2 



\ (6 - 8) 

+DJ1 + t (±x 3 - \) - \x 2 - \tu 3 + \ (xu 2 +tu 2 x ) -u x + (x- tu)u xx ) = 0, 



D t (lx 3 {3t - 1) + ±u 3 + \{u - x)u a 

+D X ( ±(u 4 - x 4 ) + h(u t (u x - 1) - (u - x)u tx + u 2 u xx + u 2 xx ) ) = 0. 



The above example illustrates that normally, to obtain simpler flux expressions, one should choose 
U = 0, unless integrals in the formula (|4.20p diverge. 

3. The scaling symmetry formula. We now reconstruct the conservation laws of the KdV 
equation (|4.16p arising from local multipliers (16, 5p using the scaling symmetry method (Section 
I4.4p with the scaling symmetry (|4.32[) . In the corresponding routine, we need to specify the point 
symmetry (|4.32p in the standard form: 

gem_get_CL_f luxes (Lam_sol , method=" Scaling" , 
symmetry={xi_x=x,xi_t=3*t ,eta_U=-2*U}) ; 
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This yields the following four conservation law expressions 

D f (2u + 3tu t + xu x ) + D x ^2u 2 + 3tuu t + xuu x + Au xx + 3tut xx + xu x 

B t {{2u + 3tu t + xu x )uj + (u(3tu txx + xu xxx ) - 3tu x u tx 
+3(2u + tut)u x x + u 2 (2u + 3tut) + u x (xu 2 — 3u x ) ) = 0, 



Df ( (2u + 3tut + xu x )(x — tu) 

+T) X ( 3tu 2 + 2u 2 (x — tu) + [xu(x — tu) — 3]u x + 3tu(x — tu)ut + 3t(tu x — l)ut x 



3{2ut + t rut — x)u xx + {x — ut)(xu xxx + 3tut xx ) = 



(6.9) 



+ 3tu t + xu x ) (\u 2 + u xx ) ) 
(j}U 3 (xu x + 3tu t ) + u 2 {3tu t + xn x - -u) + n 4 - 3tuu x u tx 
+u[xu x + 6(tn t + u)]^ + Au 2 xx + |t(n 2 + 2u xx )ut xx 
+ [^xu 2 - 3(tn tx + u x )] u xxx + (3tn t + xu x + 2u)u xxxx ) = 0, 



equivalent to the sets of conservation laws (|6.7p and (I6.9() but having significantly more complicated 
forms of fluxes and densities. 



7 Conclusion 

The direct method of construction of conservation laws is a systematic procedure applicable to all 
PDEs and PDE systems involving sufficiently smooth functions. The direct method builds on ideas 
present in the famous Noether's Theorem, in particular, the use of Euler differential operators, 
which identically annihilate any divergence expression. However, unlike Noether's Theorem, the 
direct method is not limited to PDE systems arising from some variational principle (i.e., self-adjoint 
PDE systems). 

Within the the direct method, one seeks multipliers, such that a linear combination of the PDEs 
of the system system taken with these multipliers yields a divergence expression. Importantly, for 
PDE systems that can be written in a solved form (|2.2p . the direct method is complete, since for 
such systems, all conservation laws arise as linear combinations of equations of the system. 

Determining equations for unknown multipliers are obtained from the action of Euler differential 
operators on a linear combination of the PDEs of a given system, and are in many ways similar to 
local symmetry determining equations. 

We note that the direct method can be applied without modification to ODEs and ODE systems, 
yielding their integrating factors and first integrals. 

After the multipliers are computed, fluxes of the corresponding divergence expression can be 
reconstructed using different approaches. Each approach can be better in one particular situation 
and worse or non-applicable in another (Section [5]). 

The algorithms for the direct method and the flux reconstruction have been implemented in the 
symbolic package GeM for Maple, which was used to illustrate the four methods presented in the 
present paper. 
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A Proof of Theorem [5] 

Proof. We start with a set of local conservation law multipliers {A (T [C/]}^L 1 of a given PDE system 
(j27Lj) . These multipliers satisfy (j37L]) . 
The linearizing operator associated with the PDE system ()2. 1 j) is given by 



(L R y p [u}v>> = 

a = l,...,N, 
and its adjoint by 

{h* R r P [u\w a = 



air\u] . ^[\ + ... + £M Dii ... D 



dUf dU? 



XTjP " *1 " " " % k 

h—ik 



VP, 



(A.l) 



dR a [U] 
dUP 



8R a [U] 

~dUT 



I Qp& \\J\ 

+ ... + (-l)*D il ...D ifc ( ~7^fT(r W a I , p = l,...,m, 



(A.2) 



u— ** 



acting respectively on arbitrary functions V = (V 1 ^), . . . , V m {x)) and = (Wi(x), . . . , Wn(x)). 

For each multiplier A CT [?7] = A CT (x, U, dU, . . . , d l U), introduce the corresponding linearizing oper- 
ator 



dK[U] . dK\U] , dA. a \U\ _ 



dUP dUf 



(L A ) ap [U]VP -- 
a = l,...,N, 
and its adjoint 



VP, 



(A.3) 



dUP 

+ (-l) fc D il ...D 



dUf 



(A.4) 
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W ff , p=l,...,m, 



acting respectively on arbitrary functions V = (V 1 (x), . . . , V m (x)) and W = (W\(x), . . . , Wjy(x)). 

It is straightforward to show that the operators (jA.ip - (|A.4j) satisfy the following divergence 
identities: 



W CT (L K W]^ - ^(L^)'[£7]W ff = Di&lV, W; R[U]], 



(A.5) 



^(La)^^]^ - V p (L A ) ap [U]W° = BiS l [V,W;A[U]], 

where operators S\ 5* are defined by (|4.2ip and f|4.22[) . 
To avoid the appearance of singular integrals, we need the following construction. Let 



(A.6) 



U {X) = U+(X-1)V, 



(A.7) 
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where U = (U 1 (x) , . . . , U m (x)) and V = V (x), . . . , V m (x)) are arbitrary functions, and A is a 
scalar parameter. Replacing U by Ur x \ in the conservation law identity (|3.ip . one has 

|-(A CT [[/ (A) ]iT[[/ (A) ]) = §f>i&[U {x) } = D, (J^*P(x)\) • (A.8) 

(The last identity holds since A can be viewed as an additional independent variable.) The left-hand 
side of (|A.8p can then be expressed in terms of the linearizing operators (|A.lj) and ()A.3|) as follows: 

■^(K[U W ]B7[U W ]) = K[U {x) ](L R y p [U {x) }VO + Rr[U w ](L A ) ap [U w ]V. 

Prom COj) and fOl) with W a = A a p w ] and W a = R a p( X )}, respectively, and using the 
property 

(L* R r p [U]A (7 [U] + (L* A ) ap R^[U} = 
of the adjoint operators (|A.2|) . (|A.4j) . one can show that 

j^(K[U W ]ir[U(x)]) = V(L^ p [U (x) ]K[U ( x)] +n i S i [V,A[U (x) ];R[U (x) }} 

+y(LX) CTp [t/(A)] J R' 7 [t/(A)]+D J ^[y,i?[C/ (A) ];A[C/ (A) ]] (A.9) 
= D, (s*[V,Ap {x) };Rp {x) ]} + S*[V,Rp {x) }Ap {x) }} 



Comparing (|A.8[) and (|A.9p . one finds that 

D * (^W) = D * ( si [V,HU(X)];Rp(x)}} + S i [V,R[U w ];A[U w ]] 



which implies 

^^[C/ (A) ] = 5 4 [F,A[C/ (A) ]; J R[[/ (A) ]] + 5 i [y,^[[/ (A) ];A[[/ (A) ]], (A.10) 

up to fluxes of a trivial conservation law. Now let V = U — U, for an arbitrary function U = 
(U x (x), U m {x)). Then [/ (A) = \U + (1 - X)U. Integrating ([ATTO]) with respect to A from to 
1, one obtains the desired expression (|4.20p . □ 

B Computation of fluxes of a conservation law from a symme- 
try/adjoint symmetry pair 

Suppose a given PDE system (12. ip has a local symmetry with generator (in the characteristic form) 

X = ^[u]A (B .l) 

The symmetry components fj p [u] are solutions of the symmetry determining equations, i.e., the 
linearized system 

I£[«]7H«] = 0, (B.2) 
with lT p = (Lr)' given by $J$. 
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Similarly, let {lj ct [ , u]}^L 1 be some solution of the adjoint linearized system 

L; ff [«W«] = 0, (B.3) 

with L*°" = (LJj)p given by (|A.2|) . (Functions {w ff [n]}^ =1 are often called adjoint symmetries). 
The linearizing operator (|A.1|) and its adjoint (|A.2p are related through a well-known formula [16] 

W ff L*[U]V p - V p L* p a [U}W a = Bi^ l [U} (B.4) 

where are given by 

k-ik-p-l / «Do-rrn \ 

*M = E E (°n • • • D ^) % • • • D i9 ^ 7 , (B.5) 

where ji ■ ■ ■ j q and i\. . .i v are ordered combinations of indices such that 1 < j\ < . . . < j q < i < 
i\ < ■ ■ ■ < i p < n. 
The following result is immediately obtained. 

Theorem 7. For a given PDE system (|2.ip . let be components of its local symmetry 

(jB.ip . and {w (r [it]}^ r =1 a solution of the adjoint linearized system (|B.3p . Then formulas ()4.29p with 
yp = r] p \U] and W a = oj a \U] yield fluxes of local conservation laws of the PDE system (|2.ip . 

The proof follows from noting that for the symmetry /adjoint symmetry pair {^[^j^Li, {w ff [n]}^ =1 
using V p = t) p [u} and W a = Lo a [u] in formulas (|B.4p yields Dj^ftt] = 0, which is a local conservation 
laws of the PDE system (j2"7Tj) . 

The following remarks are important. 

1. Since every set of local conservation law multipliers {A (T [?7]}^L 1 of a given PDE system (12. lh 
satisfies the adjoint linearized equations (IB.3p . it follows that in Theorem [71 one can use any 
local symmetry and any set of local conservation law multipliers of a given PDE system to 
generate its local conservation law. (Note in some cases, it is possible that a so generated 
conservation law is trivial.) 

2. A conservation law with fluxes ^[t/] (|4.29p . obtained from a symmetry/multiplier pair using 
Theorem [7J is generally inequivalent to the conservation law (|3.2p with fluxes which 
corresponds to the multipliers {A a [U]}^ =1 . One obtains the fluxes corresponding to the 
multipliers only in the case of a scaling symmetry, as discussed in Section 14.41 
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